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1 Introduction 

Summation formulas have played a very important role in analysis and number 
theory, dating back to the Poisson summation formula. The modern formulation 
of Poisson summation asserts the equality 

(1.1) E = E f/W = / fi^)e-'^"'dx) , 

valid (at least) for all Schwartz functions /. Let us take a brief historical detour 
to the beginning of the 20th century, before the notion of Schwartz function had 
been introduced. The custom then was to state (|1.1|) for more general functions 
/, such as functions of bounded variation, but supported on a finite interval, 
and usually in terms of the cosine: 



(1.2) 



E ^^^^^ ^ / ^^^^ dx + 2 E / /(^) cos(27rna;) dx ; 

a<7i<b 71=1-''^ 



the notation signifies that at points n where / has a discontinuity ~ including 
the endpoints a, b - the term /(n) is to be interpreted as the average of the left 
and right limits of f{x). Indeed, the general case of (|1.2|l can be reduced to the 
special case of a = 0, 6 = 1, which amounts to the statement that the Fourier 
series of a periodic function of bounded variation converges pointwise, to the 
average of its left and right-hand limits. 

To Voronoi and his contemporaries, formulas like H1.2|l appeared tailor-made 
to evaluate certain finite sums. That seemed significant because several impor- 
tant questions in number theory involve finite sums of arithmetic quantities. 
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In this connection, let us mention two classical examples, Dirichlet's divisor 
problem and the Gauss' circle problem. The former estimates the size of d{n), 
the number of positive divisors of an integer n. Dirichlet proved 



X 



(1.3) D{X) =dcf '^(") - ^log^-(27-l)^ = 0(^1/2), 

n=l 

where 7 is Euler's constant « .57721566. Gauss' circle problem concerns the 
average size of 

(1.4) r2{n) =dcf #{ ix,y)eZ^\x^+y^ = n}, 
for which Gauss gave the estimate 

yX 



(1.5) A{X) = 



dof 



^ — 



ttX 



Each problem has a geometric interpretation, with D{X) counting lattice 
points in the region {x,y > 0, xy < X}, and A(X) lattice points in the 
disc {x"^ + y^ < X}. These two bounds are related, as we shall see, and come 
from fairly elementary considerations. 

In the series of papers |27H29| Voronoi developed geometric and analytic 
methods to improve both Dirichlet's and Gauss' bound. Most importantly in 
retrospect, he generalized (|1.2ll by allowing weighted sums, at the expense of 
introducing more general integral operations on / than the Fourier transform. 
His formula for the divisor function reads as follows: 



, , a<n<b 

(1.6) 



d{n)f{n) = / /(.T)(logx + 27)dx + 
Y d{n) / f{x) (4i^o(47rV^) - 2^TYa{A^:^f^)) dx; 

n=l 



here / can be any piecewise continuous, piecewise monotone function, and Yq , 
Kq denote Bessel functions; for details see [T2]. Voronoi also asserted a formula 
for r2(n), which was later rigorously proved by Hardy [J] and Sierpihski I25j : 

(1.7) Y '^2(n)/(n) = Y ^2('^) / f{x)TTM2TTV^)dx, 

a<n<b n=0 

under the same hypotheses on the test function /, and with Jq referring to the 
J-Bessel function. 

These formulas have been used to improve the error estimates in (|1.3|l and 
H1.5|l to 0{X^/^) and beyond, and also in the opposite direction, to show that 
D{X) and A(X) cannot be bounded by a multiple of X^/*. Hardy and Landau 
made that observation _8,, which led them to conjecture that 0{X^/^~^'^) is the 
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best possible estimate for both the Dirichlet and Gauss problems. We shall 
return to this topic in section 

Voronoi saw H1.2|l and H1.6H1.7|I as fitting into a general pattern: he postulated 
the existence of analogous formulas for sums X]'a<n<b ^nfin), corresponding to 
any "arithmetic" sequence of coefficients a„ , 1 < n < oo . The wide latitude 
in the choice of test functions /, including characteristic functions of intervals, 
seemed like an important feature of these summation formulas. The current 
point of view is different. Characteristic functions and other discontinuous test 
functions can readily be approximated by smooth functions, and little infor- 
mation about the left hand side of (|1.2|) or 1)1. 6111. 7|l is lost by doing so. On 
the other hand, the integral transforms of the smooth approximations decay 
rapidly, which makes the right hand side of the summation formulas much eas- 
ier to analyze. We shall demonstrate this by example in section [SJ where we 
show how to obtain the estimate 0(X^/'^+'^) in Gauss' circle problem directly 
from Poisson summation. 

In an another direction, Poisson summation gave birth to modular identities, 
such as the ^-identity 

(1.8) e{t) =dcf = \e{i/t) (t>o), 

which reflects the fact that e~^^^ is its own Fourier transform. Other modu- 
lar forms, and also L-functions, can be created from Poisson's and Voronoi's 
formulas. Indeed, Riemann developed the theory of his zeta function starting 
from l|1.8() . to which he applied the Mellin transform. There arc implications 
in both directions: nowadays the coefficients of L-functions are regarded as the 
natural class of coefficients (an)n>i for which Voronoi summation formulas can 
be proved; in fact, the summation formulas are deduced from properties of the 
_L-function in question. 

The mechanism of deriving Voronoi summation from L-functions is well un- 
derstood for modular forms and Maass forms on the upper half plane, but runs 
into difficulties for more general types of automorphic forms. Recently we were 
able to overcome these difficulties by representation-theoretic arguments: our 
new result is a Voronoi formula for automorphic forms on GL(2i), which we shall 
state in section|5l It seems clear that our methods extend to automorphic forms 
on GL{n), for any n. Interestingly, the first application of our GL(3) formula 
involves Maass forms on GL{2), via the symmetric-square lift. We shall discuss 
the relevant result, due to Sarnak and Watson, following the statement of our 
formula. 

During the genesis of this note, we learned much about the history of the 
Voronoi summation formula in conversation with several colleagues. Our thanks 
go especially to Noam Elkies, Martin Huxley, Henryk Iwaniec, Philippe Michel, 
Peter Sarnak, and Trevor Wooley. 
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2 Poisson summation and the ^-function 

Andre Weil has written a fascinating historical account of the early history of the 
Riemann zeta function |31II32| . At the formal level, it turns out, the connection 
between the Poisson summation formula and the function (^(s) — X^^i 
was recognized even before Riemann. For the statement of Riemann's celebrated 
result, let us recall the definition of the related function 

(2.1) as) = 7T-^/'r{s/2)as) 

and the notion of order of an entire function /(s) - the least positive number 
fi such that |/(s)| — 0(el*l'' '), for every e > 0. 

2.2 Theorem. (Riemann |22p The function ^(s) has a meromorphic contin- 
uation to C, which obeys the functional equation ^(s) = ^(1 — s) . It has simple 
poles at s = and s = I, and the product s(s — l)'C(s) is entire, of order 1. 

In this section we shall explain the "equivalence" between theorem 12. 21 and 
the Poisson summation formula for Schwartz functions. Both results stand on 
their own, of course; the point is that each can be derived from the other. In 
one direction, Riemann proved theorem 12.21 using Poisson summation. More 
precisely, he applied the Mellin transform to the ^-relation H1.8|l . which follows 
from the Poisson summation formula with f{x) = e~^*^ ; details can be 

found in for example. 

The opposite implication is classical, too, but less well-known. We sketch the 
argument because it can serve as model for the derivation of Voronoi summation 
from properties of the corresponding L- function. First we observe that it suffices 
to deduce from theorem 12. 21 for even functions only, since both sides of the 
equation vanish when / is an odd function. In other words, we must establish 
the identity 

(2.3) "=i ° "=i 

where F{y) — i f{x) cos{2TTxy)dx , 

for every even Schwartz function /. In that case the Mellin transform of /, 

/•oo 

(2.4) Mf{s) = / fit)t^-'dt, 







is holomorphic in the region Res > — 2 , except for a simple pole at s = 0, with 
residue /(O). The Mellin inversion formula 

(2.5) fix) ^ ^1 Mf{s)\x\'^ds (a>0) 
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allows us to express the left hand side of (|2.3|) as 

oo 

(2.6) 5^ /(n) = — / as)Mfis)ds (a > l) . 

n=l "'R<= s='^ 

When we shift the contour of integration in this formula from Re s = a > 1 to 
Re s = — 1 , we pick up a pole at s = 1 with residue M/(l), and a pole at s = 
with residue C(0)/(0) — —1/2/(0). These are the first two terms on the right 
hand side of H2.3|l . Note that the contour shift is justified: ^(s) grows slowly 
and Mf{s) decays rapidly as |Im s| ^ cx) . We now change variables from s to 
1 — s and apply the functional equation: 

f{n) M/(l) + M f ^(,) Af/(,s) ds 

(2.7) 



2m Jncs=2 n('^-^)/^rC-^) 



1 f TT-'^^TiD 
1 



^ C(s) M/(l -s)ds. 



But 

(2.8) MF{s) = ___l^M/(l-s) 



is the Mellin transform of the function F in the identity H2.3|l . This follows from 
the classical identity 

(2.9) / cos(2^x) \xr' dx = ,^_^.2la-s, (0 < Res < 1), 

initially in the region {0 < Res < 1}, but then for all s £ C by meromorphic 
continuation. The even Schwartz function F can play the part of / in (|2.()|l . At 
this point, (|2.7I - I2.8|) imply ifT^ . 



3 L-functions and Voronoi summation 

The two Dirichlet series with coefficients d{n) and r2 (n) , respectively, also sat- 
isfy functional equations. In the case of d(n), 

(3.1) 5]°° d(n)n-^ = ^A^^r' = C'{s) 

^ — ^n— 1 ^ — ^m.n>l 

inherits its functional equation and analytic properties from the Riemann zeta 
function. The generating series for r2(n), 

(3-2) Cq(v^)(s) = jEr=i'^2(")""'' 
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is known as the Dedekind zeta function of the number field 1), because it 

counts the number of ideals of norm n. It factors as the product of two Dirichlet 
series: 

(3-3) Cq(V^)(s) = Cis)L{s,Xi), 

where 

(-1)'= 



fc=o {2k + ly 

is the first example of a Dirichlet L-fimction. It is entire, of order 1, and satisfies 
the functional equation 



(3.5) L(.,X4) = 2^-" ^_(.^i)/,r(4i) Lil-s,X4). 



The factorization H3.3|l is a result from class field theory, which reflects the fact 
that an odd prime can be expressed as the sum of two squares if and only if it 
is congruent to 1 modulo 4. The factorization is also equivalent to the identity 

r^in) = 4 e N I d|7i, d=lmod4} 
''^'^^ - 4 #{de N I d|n, d= -lmod4}, 

which implies the bound 

(3.7) r2(n) < 4d(n) = 0{n') , for any e > ; 

this uses a standard estimate for d{n), which can be found in |^ p. 186], for 
example. 

The Voronoi summation formulas H1.6|l and (|1.7|) can be derived from the 
functional equations of the Dirichlet series H3.1|l and H3.2|l . much in the same 
way as Poisson summation was derived from theorem 12.21 To match the F- 
factors in the functional equations to the integral kernels on the right hand 
sides of p.6H1.7|l . one uses the identities 

{4Ko{Vx)~2TTYQ{y^))x-'^^dx = 



(3.8) = 2 



Jo{Vx)x^ ^dx 



4s 2s 
TT 



m-s) 



which are valid for < Res < 3/4. In effect, these identities take the place of 
(|2.9|l in the argument in the previous section. 

Modular forms are a rich source of L- functions, and each such L- function 
satisfies a functional equation from which a Voronoi formula can be deduced. 
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To simplify the discussion, we shall consider the case of "full level" , i.e., modular 
forms invariant under all of S'L(2,Z). Recall the definition of modular form of 
weight k: a holomorphic function F : {Imz > 0} — > C such that 

(3,9) ^w^'^+^'^KS^) 'f 

and \F{z)\ = 0((lmz)^) as Imz^oo, 

for some iV > 0. The invariance condition implies invariance under z i-^ z + 1 in 
particular, so F{z) has a Laurent expansion in powers of e'^^*^. The condition of 
moderate growth rules out negative powers. We may and shall exclude even the 
0-th power - in other words, we shall assume F is cuspidal. Eisenstein series, 
which span a linear complement of the space of cuspidal modular forms, require 
a separate but simpler analysis. In the cuspidal case, then. 



Eoo 
n— 1 



^(fc-l)/2 g27rinz 



Traditionally one denotes e^'^*^ by the symbol q and calls H3.9() the q-expansion, 
but we shall not follow that practice. 

The Dirichlet series with coefficients a„ is called the standard L-function of 
the modular form F, 

Eoo 
n—l 

It satisfies a functional equation, which reflects the transformation law for F 
with respect to z i-^ — 1/z. We should remark that the factor Inl'*^^^)/^ in (|3.1U|) 
has the effect of making the functional equation relate L{s, F) to L{1 — s, F) 
rather than to L{k — s,F), as is the case in the classical literature. The L- 
function (|3.11l) can be "twisted" , either by an additive character x of finite 
order, 

(3.12) x{n) = e^""'""^", with a,ceZ, c 7^ , (a, c) = 1 , 
or by a Dirichlet character modulo q, 

(^) ^ / Xq(N) if {n,q) = 1 

(3.13) ' \ a {n,q)^l, 

for some character Xg ■ (Z/^/Z)* C* , 
with [n] denoting the image of n in (Z/qZ)*. The Dirichlet series 

(3.14) L^{s,F) = J2 1 Xin)a^n-^ 

^ — ^n— 1 

is the twisted L-function of F, with "additive twist" in the case of H3.12|l . 
respectively "multiplicative twist" in the case of H3.13|l . Multiplicative twists 
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have arithmetic significance. For example, multiplicatively twisted i-functions 
of Hecke eigenforms have Euler products and are considered special cases of 
Langlands L-functions. Additive twists, on the other hand, play an important 
role in the analytic study of i-functions. 

A Dirichlet character x raodulo q is called primitive if, in the notation of 
Xq does not arise as the puUback of a character of (Z/g'Z)* for some divi- 
sor q' of q. In this situation, the multiplicatively twisted L-function (s, F) also 
satisfies a functional equation. Everything that has been said so far can be gen- 
eralized in some fashion to automorphic forms on GL{n) and other higher rank 
groups. Not so the assertion that the additively twisted L-function L^{s,F) 
satisfies a functional equation, which is true for GL{2) but fails already for 
GL(3). 

Maass forms are eigenfunctions of the Laplace operator on the upper half 
plane, of moderate growth, invariant under the natural action of SL{2,'Z) or 
more generally, invariant under a subgroup F C 5'L(2,Z) of finite index. Just 
as in the case of modular forms, one can attach L functions to Maass forms. 
These L-functions, both with multiplicative and additive twists, again satisfy 
functional equations. 

In 1930, Wilton proved a Voronoi summation formula for the modular form 
A, the essentially unique cusp form of weight 12. Wilton's argument applies to 
any cuspidal modular form, even with additive twist. It can also be adapted to 
Maass forms. The reason, of course, is the functional equation for the additively 
twisted L-functions, which is used in the same way as the functional equations 
for C(s), C^('S), and Cq(^/^) (s). This type of argument does not extend to groups 
beyond GL{2), however, since then additively twisted L-functions no longer 
satisfy functional equations. Recently we were able to prove a Voronoi formula 
for GL(3) with representation-theoretic arguments flM- It seems clear that our 
argument works for GL{n). 

In the next section we state the Voronoi summation formula for modular 
forms and Maass forms and then sketch a proof. Our argument is essentially 
different from the usual arguments, in that it works directly with automorphic 
forms; in effect, it is an adaptation to GL{2) of our argument for GL(3). We go 
on to discuss our GL(3)-formula in section|Sl 



4 Voronoi Summation for GL{2) 

Cuspidal Maass forms and cuspidal modular forms for F = SL{2, Z) correspond 
to embeddings of, respectively, principal series and discrete series representa- 
tions into L^(F\S'L(2, M)). Let us begin with some quite general remarks about 
discrete summands in L^(F\G). 

Initially G will denote an arbitrary unimodular Lie group of type I, F C G a 
discrete subgroup, and to : Zq ^ C* a unitary character of Zq = center of G. 
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Then G acts unitarily on 

Ll{T\G) = { / G LL(r\G) I Jigz) - u{z) f{g) 

for ze Zg, and Jy\g/Zg '-^'^ c^ff < oo } . 

If (tt, y) is an irreducible unitary representation of G and 

(4.2) j:V ^ Ll{r\G) 

a G-invariant embedding, the subspace V°° C F of C°° vectors gets mapped 
to F-invariant smooth functions on G, which can be evaluated at the identity. 
Thus 

(4.3) T, ^ C, T,{v) = jiv){e), 

is a well defined, F-invariant linear map. Since j(u)(e) can be bounded in terms 
of finitely many derivatives, one can show that Tj is continuous with respect 
to the natural topology on . In short, tj can be regarded as a F-invariant 
distribution vector for the unitary representation (tt', V') dual to (vr, V), 

(4.4) r, e(y'-^)^ 

We shall call Tj the automorphic distribution corresponding to the embedding 
j. It determines j completely since j{v){g) — j(7r(g)w)(e) = Tj{Tr{g)v) for any 
V in the dense subspace V°° C V. Our discussion will focus on Tj. To simplify 
the notation, we switch the roles of (tt, V) and (tt', V'); from now on, 

(4.5) j-.V ^ ^^i(r\G), r, e(y— )r. 

This is legitimate since integration over T\G / Zq sets up a nondegenerate pair- 
ing between Ll{T\G) and Ll_,{T\G). 

We now specialize our hypotheses to the case of G = S'L(2,M) and F = 
S'L(2, Z). Since F contains Zq: only discrete series representations with even 
minimal weight and even principal series representations can occur as sum- 
mands in (|4.4|l . Discrete series representations have realizations as subrepresen- 
tations of non- unitary principal series representations. According to a result of 
Casselman and Wallach - whose specialization to G = SL(2, K) can be ver- 
ified directly - the assignment V i-^ V~°° is an exact functor. Thus, even when 
TT is a representation of the discrete series, the automorphic distribution Tj may 
be regarded as a distribution vector for a (possibly reducible) even principal 
series representation. As model for even principal series representations we take 

(4.6) = {f eC^{R) I Ixl^-'-'fil/x) eC°°{R)}, 
with J/ G C as representation parameter. The group G acts via 

(4.7) (nAg)f){x) ^Icx + d^-'fi^J if 
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The resulting representation is irreducible whenever ^ Z + 1/2 ; in that case 
TT^ ~ TT-i, .If = 1/2 — fcgZ + 1/2 is negative, contains two discrete series 
representations as subrepresentations and has a {2k — l)-dimensional quotient, 
whereas has a {2k — l)-dimensional subrepresentation and two discrete 
series quotients. 

In contrast to H4.5|l . the restriction map from the space of distribution vectors 
V~°° to C~°°(M) is not injective: a distribution on a manifold is not determined 
by its restriction to a dense open subset. Still, in terms of the restriction map 

(4.8) V-°° — > C-°^\R), 

the action of G is given by the formula H4.7|) at all points a; G R such that 
cx + d^O. 

We now suppose that tj e (V'j7°°)'" is a cuspidal automorphic distribution - 
i.e., an automorphic distribution corresponding to a cuspidal embedding H4.5|) . 
The automorphy with respect to F implies in particular that Tj is periodic of 
period 1. Consequently its restriction to R has a Fourier expansion 

(4.9) T,{x) = ^ ajnr^'e^-'-, 

whose constant term vanishes because of the assumption of cuspidaUty. In the 
case of an embedding of a discrete series representation, a„ — for either all 
positive or all negative indices n. In the case of a principal series representa- 
tion, the parameter v is determined only up to sign. The two possible choices 
determine different distributions in C^°°(R), which are related by the standard 
intertwining operator. One can show that the coefficients a„ corresponding to 
the two choices agree up to a multiplicative constant - that is accomplished by 
the factor \n\~'' in H4.9(l . If tj satisfies the conditions of automorphy under F 
and cuspidality, then so does x i-^ Tj{—x). For discrete series representations, 
this operation switches the roles of the two subrepresentations in the ambient 
V~°° . In the Maass case, on the other hand, this operation allows us to con- 
sider even and odd automorphic distributions separately. In other words, in the 
Maass case we may and shall assume 

(4.10) Tj{-x) = {-l)'^Tj{x) {r]eZ/2Z). 

When this condition is imposed, the a„ for n > determine also the a„ with 
negative index, and 

Eoo 
n— 1 

is the standard L-function of the cuspidal modular form or Maass form corre- 
sponding to the embedding H4.5|l : one sees this most easily by identifying the 
action of the Hecke operators on the Fourier series 1)4. 9|) . 

For the statement of the Voronoi formula for GL{2), we fix an automorphic 
distribution tj as above, as well as integers a, a, c such that c 7^ 0, (a, c) = 1, 
and da=l modulo c. 
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4.12 Theorem. pi llllll5| If the Schwartz function f E S{M.) vanishes to infi- 
nite order at the origin, 



Care must be taken in interpreting the integral expression for F. When the 
integration is performed in the indicated order - first with respect to xi, then 
with respect to X2 - the inner integral converges absolutely, but the outer in- 
tegral only conditionally, and only when Rez/ > 0. For other values of v, the 
integral retains meaning by analytic continuation. Alternatively F can be char- 
acterized by a formula for its Mellin transform in terms of the Mellin transform 
of /; this is the form in which the theorem is usually stated. The proof shows 
that F and all of its derivatives decay rapidly as \x\ — > oo. At the origin F has 
certain mild singularities; see OI for details. 

The Voronoi summation formula for GL(2) has become an essential tool in 
analytic number theory, for problems such as sub-convexity of automorphic L- 
functions or the analysis of the Petersson and Kuznetsov trace formulas. In all 
the applications the presence of the additive twist plays a critical role. A formal 
application of the theorem with f{x) = (sgna;)''|a;|~* - which is not covered by 
the hypotheses, of course - gives the functional equation for the i-function with 
additive twists. In fact, the proof that we shall sketch can be adapted to deal 
even with this "illegitimate" choice of test function. The functional equation 
for multiplicative twists can be deduced by taking appropriate linear combina- 
tions of the additively twisted functional equation. Conversely, the functional 
equation with multiplicative twist implies the equation with additive twist, and 
that in turn can be used to establish the summation formula, along the lines of 
the arguments in section 1^ pi lTTlIT^] . 

For GL(3), the passage from multiplicative to additive twists is far more 
subtle; in fact, additively twisted L-functions for GL(3) do not satisfy functional 
equations, as was pointed out already. This has been the obstacle to a Voronoi 
formula for GL(3) in the past. Our proof of the GL(3) formula in overcomes 
the obstacle by working directly with the automorphic distribution. We shall 
now sketch a proof of theorem 14.121 by adapting our argument in TU" to the 
case of 5'L(2,M). 

Because of the hypotheses of the theorem, there exits b G Z such that aa — 
1 -|- be. The 2x2 matrix 7 with entries a, 6, c, a lies in T. We now use the 
formula H4.7I) to make the automorphy of Tj under 7 explicit, then substitute 
X — a/c for X, which results in the equation 




with F described by the conditionally convergent integral 




dxi dX2 . 



(4.13) 
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We integrate this equality of tempered distributions against a test function g 
in the Schwartz space 5(R). On one side, using the Fourier expansion 1)4. 9|l . we 
get 



(4.14) 



'^3 (a; - f ) 9(x) dx = 

= E ^ a„|nr''e"2-^"'^/^?(-n). 



The other side of H4.13|l appears to be ill defined at a; = 0. However, the fact 
that Tj is periodic without constant term implies that this distribution can 
be expressed as the iV-th derivative of a bounded continuous function for ev- 
ery sufficiently large integer iV. This fact, coupled with an integration-by-parts 
argument, can be used to show that the distribution Tj(1/x) has a canonical 
extension across the origin. With this interpretation, (|4.13|) remains valid even 
at a; = 0; roughly speaking, the cuspidality hypothesis rules out the delta func- 
tion or derivatives of the delta function from contributing to the expression for 
of Tj at infinity. In short, the formal computation 

ca;p--i r, (f - ^) g{x)dx = 

V kP''^^ ar.W]-" e^'"'"'"^'^'"'-^'''^ g{x)dx 



(4.15) " 1^1'" 



can be made rigorous |2Uj . 

Let / e 5(R) be a Schwartz function which vanishes to infinite order at x = 0, 
as in the statement of the theorem. Then g(x) = J^^f {y)\y\'^ e{—xy)dy is also a 
Schwartz function, and f{x) = \x\~" g{—x). With this choice of g, (|4.13ll4.l3|l 
imply 



n#0 



1^12.-1 



^ JJr2 
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Setting t — nc ^ and changing the symbols for the variables of integration we 
get the statement of the theorem. 



5 The GL{3) summation formula 

The discussion in the beginning of the previous section applies also to the 
setting of G = GL(3,R), F = Gi(3,Z). Now G has one dimensional center, 
so the central character uj in (|4.1H4.5|) is necessary - without it, the space of 
F-invariant functions would vanish. 

As in the case of S'L(2,R), every cuspidal automorphic distribution can be 
regarded as a distribution vector for some principal series representation; this 
assertion depends on the Casselman embedding theorem and the Casselman- 
Wallach exactness theorem for spaces of distribution vectors. Not-necessarily- 
unitary principal series representations of G with unitary central character cor- 
respond to parameters (A, (5), 

A = (Ai, A2, As) e C^, 5 = (<5i, (52, ^3) e (Z/2Z)3, 
E,ReA,=0. 

The Sj specify the inducing character on the component group of the diagonal 
Cartan subgroup, and the Xj are normalized so that the Weyl group action fixes 
A = and p = (1,0,-1) is the half-sum of the positive roots. The condition 
Re Xj = reflects the requirement that Zq must act unitarily. As model for 
V^g we take the space of smooth functions on 













(5.2) 


7V+ = < 


1 B 


X, y, z £ M. 


> C GL(3, 













with an appropriate regularity condition at "infinity" and an action ttx^s of 
G involving a factor of automorphy. In analogy to l|4.8|) . there is a natural 
restriction map 

(5.3) i/,:r c-°^{N+). 

The subgroup C G acts by translation, without factor of automorphy, via 
this restriction map. 

We now consider a particular cuspidal automorphic distribution tj G ^xs"- 
By restriction to N^, and in terms of the coordinates x, y, z in H5.2() . we may 
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think of Tj as a distribution in these variables. Naively speaking, the restriction 
to does not determine Tj, but the F-automorphy of Tj ensures that no 
information is lost. The invariance under 7V_|_ H T allows us to regard Tj{x, y, z) 
as a distribution on the compact manifold {N^ n r)\A^+. As such, it has a 
Fourier decomposition with respect to the right action of A^+t 

(5-4) _ V- |n|^-^(sgnn)^^ ^ 2^^i„^x+ny) , 

^ 771^1 (sgnm)*i 

here we have written out the abelian part of the decomposition explicitly, but 
the remaining components - those on which the center of Af_|_ acts according to 
a non-trivial character - are summarily denoted by . . . . The cuspidality of Tj 
accounts for the absence of abelian terms indexed by (0, n) or (m, 0). The powers 
of |m|, \n\ and the signs play the same role as the factor \n\~'^ in (|4.9|) - they 
make the coefRcients independent of the choice of Casselman embedding, up to 
universal multiplicative factors. With this normalization of the coefficients, 

oo oo 

(5.5) L{s,Tj) = ^ai,„77"'', L{s,Tj) = a^^i 

n—l n—1 

are the standard L-functions of the automorphic representation corresponding 
to Tj and the dual automorphic representation; the latter corresponds to the 
dual automorphic distribution fj whose abelian Fourier coefficients agree with 
those of Tj except for a transposition of the indices. 

The statement of our Voronoi formula for GL(i) is entirely analogous to that 
for GL{2). It involves the Kloosterman sum 

(5.6) S{m,n;c) = ^ e2'^^("^+"*)/'= {m,n,c e Z , c ^ 0) ; 

a;e(Z/cZ)* 

in the summation, x denotes an integer which represents the multiplicative 
inverse of x modulo c. 

5.7 Theorem. Let f G 5(M) be a Schwartz function which vanishes to infi- 
nite order at the origin, or more generally, a function on R — {0} such that 
(sgna;)*^|x|-^3/(a;) g ^(r). Then for (a,c) = 1, c ^ 0, aa = I (mode) and 

q>0, 



n^o \n\ \ c?q 



i\cq 

where 



Q 1 — Ai — A2 — A 

d 



F{t)^ I /(^i^)n(e-2--.>,.|-^.(sgna;,)^^)dx3dx2dxi. 

JR3 ^ * ^ fj^ 
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The integral defining F converges when performed as repeated integral in the 
indicated order ~ i.e., with x^, first, then X2, then xi - and provided ReAi > 
Re A2 > Re A3 ; it has meaning for arbitrary values of Ai , A2 , A3 by analytic 
continuation. As in the case of GL{2), the function F and all its derivatives 
decay rapidly as |a;| — > 00, but F has a mild singularity at the origin. One can 
also characterize F in terms of its Mellin transform |19) . 

While the statement of the theorem only involves the abelian Fourier coeffi- 
cients, the proof uses the non-abelian coefficients in an essential way. The invari- 
ance of Tj under the two simple Weyl reflections - which have representatives in 
r - translates into relations between the abelian and non-abelian coefficients. 
Additive twists show up transparently in these relations. A careful analysis then 
leads to the summation formula. Variants of these arguments can be used to 
prove the functional equations for the standard L-functions, with multiplicative 
twists, and also to prove the "converse theorem" of Jacquet, Piatetski-Shapiro 
and Shalika jl4| : an automorphic representation can be reconstructed from the 
two standard L-functions, provided they have Euler products of the appropri- 
ate type, their multiplicative twists satisfy functional equations, and they grow 
slowly in the imaginary directions jl9| . 

Cuspidal GL (3, Z)-automorphic representations exist abundantly, as a con- 
sequence of variants of the Selberg trace formula jl8| . Explicit construction 
of such representations is a very different matter - for that, one relies mainly 
on the Gelbart- Jacquet symmetric square lift of cuspidal 5*272 (Z)-automorphic 
representations Not surprisingly, then, our GL{3) formula has been applied 
first to a problem about modular forms on GL(2). Indeed, when Peter Sarnak 
saw our analysis of automorphic distributions on GL(3,M), he suggested we 
should be able to establish a Voronoi summation formula for GL(3), a formula 
he considered potentially useful for the solution of a certain problem about 
Maass forms. Let us describe this problem. 

Let X be a compact Riemann surface, {(/'j j^g an orthonormal basis of 
Laplace eigenf unctions, and fxj the eigenvalue corresponding to the j-th eigen- 
function. Sogge |2E| bounds the L^'-norms of the (f)j in terms of the fij. His 
bound comes from differential-geometric estimates and is sharp when applied 
to a general compact Riemann surface. A stronger bound is expected to hold 
for the discrete spectrum of a complete hyperbolic surface of finite volume, 
even when the surface fails to be compact. In the special case of the surface 
X = SL2{'£)\'H, Ti. = {Rez > 0}, Berry's random wave model and experi- 
mental data of Hejhal-Rackner [S] suggest that the distribution of the values of 
eigenfunctions tend towards a Gaussian distribution. In particular this would 
imply the estimate 

(5.8) UjWj, = 0{ii)) (2<p<cx3, e>0) , 

for any orthonormal basis of Hecke-Laplace eigenfunctions of the discrete spec- 
trum of the modular surface SL2{1j)\H. The survey explains these conjec- 
tures as eigenfunction analogues of the Lindelof conjecture; in fact, a variant of 
H5.8|l would imply the classical Lindelof conjecture. 
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Sarnak and Watson [21] have recently announced (|5.8|) for the case p = A - at 
present under the assumption of the Ramanujan conjectures for Maass forms on 
GL{2). In contrast, the bound of [25) for p = A, were it to apply in this noncom- 
pact situation, only predicts the bound 0{^^J^^). Sarnak- Watson's argument 
uses our Voronoi summation formula for GL(3), among other ingredients. 

To show how our formula can be applied to modular forms and Maass forms, 
we shall describe the abelian Fourier coefficients of the symmetric square lifting 
terms of the coefficients of the original modular or Maass form. Thus we consider 
a cuspidal 5'L(2,Z)-automorphic distribution r, with representation parameter 
I/, which we express as in (|4.9|l . Recall that a„ = for n < and v ~ 1/2 — fcif 
T corresponds to a holomorphic cusp form of weight 2k. In the Maass case, we 
impose the parity condition (|4.1UI) . which is equivalent to 



(5.9) 



= (-1) 



{t] e Z/2Z) . 



Still in the Maass case, the parameter v is purely imaginary, but only determined 
up to a sign factor. In both cases, the a„ with n > completely determine r. 

The symmetric square lift is defined for Hecke eigenforms. The Hecke property 
implies in particular that oi 0, making it possible to rescale r so that ai — 1. 
In that case, the Hecke property is equivalent to the Euler product factorization 



(5.10) 



L(s,r) = £, 



-s\-l/ 



o^pP ") ^(1 - "p P 
with ttp + 



l„-s^-l 



The Gelbart- Jacquet symmetric square lift of r is a cuspidal GL{3, Z)-automor- 
phic distribution Sym^ r, whose abelian Fourier coefficients we denote by am,n- 
Recall the definition of the Mobius function /i : N ^ {±1}, 



(5.11) A.(prpr---p^) 



(-1)'' if aU ej = 1 and 

Pi ^ p-j for ij^j, 







if at least one e, > 1 , 



5.12 Proposition. The lifted automorphic distribution Sym^ r has represen- 
tation parameters 



A = (2^, -2i^, 0) , 6 = 
and abelian Fourier coefficients 



(0, 0, 0) in the Maass case 
(1,0,1) in the holomorphic case 



'im.n = ^ Kd) ^ {OrOs \r 
d\{m,n) 



fc,£e N}, 
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With these choices of A, S and a„i^n, theorem 15 . 71 provides information about 
the an ■ In the Maass case, the representation nx^s of GL{3, R) which corresponds 
to the parameters (A, S) is an irreducible principle series representation, and in 
the holomorphic case, t:x,s is unitarily and irreducibly induced from a discrete 
series representation of GL{2,M.). 

Proof. The completed i-function of the automorphic representation determined 
by T is the product 

(5.13) A(s,r) = ioo(s,r) ifi„ito(s,T) , 

with £finitc('S, t) = L(s, r), as in (|5.1U|) . corresponding to the finite places and 

Loo{s,t) = 

(5.14) ^ j r^{s + + 1]) r^is - V + T]) in the Maass case 

I rc(s + fc — i) in the holomorphic case 

corresponding to the archimedean place. Here we are following the notational 
convention 

(5.15) Tr{s) = ^-^/2r(s/2), Tc{s) = 2(2^)-^r(s), 

which is commonly used in the context of L-functions. The completed i-func- 
tion H5.13|l satisfies the functional equation 

(5.16) A(s,t) = wA(1-s,t), withweC, H = I ■ 

Gelbart- Jacquet ^ characterize the symmetric square lifting Sym^ r by its 
completed L-function 

(5.17) A{s,Sym^T) = Loo(s, Sym^ t) Lfi„itc(s, Sym^ r) , 

which satisfies the same type of functional equation as A(s, r). The archimedean 
part is 

-^oo(s,Sym^T) = 

(5.18) _ J rR(s + 2i/)rR(s)rR(s - 2i/) in the Maass case 

I Tcis + 2k — l)rR(s + 1) in the holomorphic case 

and the finite part 

-^-finitc(s, Sym r) = > AnU^'^ ^ 

^ — ^n—l 

p 

Expanding the Euler product one finds 



(5.20) Lfi„ite(s, Sym^ r) = C(2s) V 

^ — 



oo 

arj2 n 
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which in turn impUes 

(5.21) An - II^.|„ «(n/<i^)^ ■ 

In |19[ Theorem 6.15] we relate the coefficients a„.„i and parameters A, 5 of a 
cuspidal Gi(3, Z)-automorphic representation of GL(3,]R) to its standard L- 
function. In the current setting this gives 

(5.22) ai,n = a„,i = yl|„| (n^^O). 

The abehan Fourier coefficients am,n of any Hecke eigendistribution for GL(3, Z) 
satisfy the product relation 

(5.23) ara,n = V"^,, , M(c^) « 1 . n/d « m/d . 1 • 

In the situation at hand this leads to the formula for the a^.n asserted by the 
proposition. The values of the parameters A, (5 in the statement of the proposi- 
tion match l|5.18|l to the calculation of the F-factors in Theorem 6.15]. 

6 Error term in the circle problem 

In this section we try to give a glimpse of how the Voronoi summation formula 
is used in practice. Voronoi proved the bound 

(6.1) A(X) = 0(Xi/3) 

for the error term H1.5|l . We shall show how to obtain the slightly weaker bound 

(6.2) A(X) = 0(Xi/3+<=) 

from fairly simple estimates. Voronoi's formula H1.7|l for the sequence r2(n) can 
be established from the Poisson summation formula in two variables, applied to 
radially symmetric test functions /. In the interest of brevity we shall deduce 
H1.7|l directly from Poisson summation in two variables, instead of taking a 
detour via (|1.7|) - which would not fundamentally change the argument. 

Let xu denote the characteristic function of I?, the unit disc in M^. In the 
following, X > \ will be a large parameter and 8 = S{X) = X^^^^. We choose 
a radially symmetric, smooth function <& > 0, supported on the unit disc D, of 
total integral 1. The convolution product 

(6.3) Fxiv) ^ S-^xd{X-^^^v)*<P{S-'^v) 

is also radially symmetric and smooth, < Fx < 1, Fx{v) = 1 for < ^/X—S 
and Fx{v) = for \\v\\ > ^/X + S. Its Fourier transform 



(6.4) 



Fx{v) = X xhiX'^^v)^Sv) 
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decays rapidly and is again radially symmetric. Because of the radial symmetry, 
there exist functions fx, tp '■ R>o ^ R such that 

(6.5) Fx{v) = fx{\\v\\), Hv) = i^iWvW). 

The Fourier transform of xd is expressible in terms of a Bessel function: 

(6.6) xhiv) = \\v\\-' M27r\\v\\) 

IHl p. 962]. Since xd and $ have total integral tt and 1, respectively, 

(6.7) xh{0) = V(0) = 1; 

the former also follows from the known properties of Ji, of course. Poisson 
summation in identifies the sum X^ugz^ ^x{v) with the corresponding sum 
for the Fourier transform. In view of (|6.4H6.7| ). we conclude 

J2 ■^^r2{n)fx{V^) ^ nX + 

^ — ^n>0 

+ V n-^^^r2{n)iP{Sy/^)Ji(2nV^). 



(6.8) 



Recall that the difference between fx and the characteristic function of the 
interval [0, Vx] is supported on the interval [\/X — 6, ^/X + S] and has absolute 
value < 1. Thus 



(6.9) 
where 



(6.10) 



A(X) < Bi{X) + B2{X) + BsiX), 

Bi{X) = ^{r2{n)\VX-5<^/^<^+5}, 
B2{X) ^ VX y i^{5^) .h{2^V^) 

l<Ji<(S-2 ^ 

r2{n) 

<n<oo 



B:i{X) = VX Y i;{5V^) Ji{2nV^) 



To establish ()6.2|l . it suffices to bound each of these three terms by a multiple 
of for every e > 0. 

In the case of Bi{X), this follows - with room to spare - from the definition 
S = g^^^ ^j^g bound r2{n) = 0{n'); cf. For B2{X), we also use the 

stationary phase estimate 

(6.11) |Ji(a;)| = 0(a;-i/2) 

IHl p. 972] and the boundedness of ip: 



(6.12) 



B2iX) ^0[X'/^J2, 



-3/4^ 



<n<5-2 
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We argue similarly in the case of B3{X), but now use the rapid decay of tp. In 
particular TpiS^) = ©(ni/^X-i/^)- V2-3e^ j^g^^^g 

(6.13) B^iX) = o(x'/'+'/^y n-'-'/A = 0(^1/3+^). 

That completes the argument. 

With more effort, one can remove the e from these bounds and get Voronoi's 
result H6.1|l . A more serious limitation of the argument above is the use of 
(|6.11ll to estimate the size of Ji. Various authors have done better by taking 
advantage of the oscillation of Ji |7ll51 [T^[TCl[r7| ; the current record bound for 
A(X) is 0(X23/73+e) [XH]. In the other direction, Hardy [7| used the oscillation 
to show that A(X) is as big as a constant times X^/'* infinitely often. That 
makes 0{X^/'^^'-) the best possible estimate for the error term, conjectured by 
Hardy and Landau [Hj in 1924. This is analogous to the Riemann hypothesis 
and the prime number theorem, where the optimal bound on the error term 
7r(X) — Li(X) is thought to be 0(|X|^/^ logX), roughly the square root of the 
trivial bound 0{X/ \ogX). However, the issues are different: the circle problem 
is one in analysis, in that it asks for a sharp bound for a sum when a harsh 
cutoff function is introduced; formula (|6.8I) already gives precise estimates for 
smoothed sums. In the case of the Riemann hypothesis, the optimal bounds are 
not known even for smoothed sums. 
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